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Abstract. Any characterization of a single-photon source is not complete without 
specifying its second-order degree of coherence, i.e., its g'^'^^ function. An accurate 
measurement of such coherence functions commonly requires high-precision single- 
photon detectors, in whose absence, only time-averaged measurements are possible. 
It is not clear, however, how the resulting time- averaged quantities can be used 
to properly characterize the source. In this paper, we investigate this issue for a 
heralded source of single photons that relies on continuous-wave parametric down- 
conversion. By accounting for major shortcomings of the source and the detectors — i.e., 
the multiple-photon emissions of the source, the time resolution of photodetectors, and 
our chosen width of coincidence window — our theory enables us to infer the true source 
properties from imperfect measurements. Our theoretical results are corroborated by 
an experimental demonstration using a PPKTP crystal pumped by a blue laser, that 
results in a single-photon generation rate about 1.2 millions per second per milliwatt 
of pump power. This work takes an important step toward the standardization of such 
heralded single-photon sources. 
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1. Introduction 

Single-photon sources (SPSs) are important elements in quantum communication, 
optical quantum computing, and metrology [Ij. To satisfy the requirements for such 
applications, it is desired that such SPSs only create single, and not multiple, photons 
in an on-demand way. To build such a source, one can employ an array of heralded 
single-photon sources (HSPS) [2], with possibly low individual probability of single- 
photon generation, but with an overall probability that approaches one for a sufficiently 
large number of sources in the array [3j. One of the most convenient ways to generate 
heralded single photons is based on spontaneous parametric down-conversion (SPDC) 
|4j. In this scheme, by pumping a crystal with x^^^ nonlinearity, one can generate 
a pair of signal and idler beams whose photon numbers in given time intervals are 
highly correlated, thus by detecting a single photon on the idler beam one can, ideally, 
guarantee the presence of a single photon on the signal beam. Inevitable to this scheme 
is the occasional generation of multiple-photon packets in each beam, which degrade 
the reliability of the SPS. This effect can be best examined quantitatively by evaluating 
the degree of second-order coherence (termed coherence function, hereafter, for brevity) 
for our source. In this paper, we analytically calculate the coherence function for our 
SPDC-based source and measure it in an experiment. We consider the impact of finite 
time resolutions on the measurement results, and how this shortcoming may affect the 
proper characterization of such sources. In fact, our work paves the way for developing 
standard specifications for HSPSs. 

For an ideal SPS, we expect that its second-order degree of coherence, viz. its 
^'(^^ function, is zero at the origin [5j. This is equivalent to having no coincidence 
detection on the two detectors of a Hanbury-Brown and Twiss (HBT) interferometer 
[6]. "Coincidence," in theory, refers to two simultaneous events. In practice, however, 
a coincidence event can only be verified within a certain accuracy permitted by the 
employed measurement devices, e.g. photodetectors and their respective electronics. 
This requires us to define coincidence by referring to two events that occur within a 
coincidence window whose width is greater than zero. The measured values for g^'^'^ — 
or even its simplified form commonly expressed as the probability of having two or 
more photons in a given (short) time interval over the probability square of having 
only one photon — may well depend on our choice of coincidence window as well as 
on other experimental parameters. Such dependence poses a challenge on the proper 
standardization of HSPSs because a single value of ^'^^^(0) does not necessarily convey 
sufficient information to characterize such a source. This is particularly the case in our 
continuous- wave (cw) SPDC-based HSPS, whose coherence functions may have widths 
in the sub-picosecond regime, much lower than what typical photodetectors can measure. 

The measurement of coherence functions in SPDC-based HSPSs is not only affected 
by the above time parameters but also by the multiple-photon contribution to the 
SPDC output. The latter is a function of the pump power, which, at the same time, 
determines the rate of single-photon generation of our source. In this paper, we present 
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a theoretical framework that not only, for the first time, accounts for the multiple- 
photon emission in our source but also allows us to examine the effect of imperfect 
devices on the g^'^'^ measurement. Such an analysis provides prescriptions for proper 
characterization of coherence properties of HSPSs and how such figures can be measured 
in practical experimental setups. This is of crucial importance because such devices have 
already been introduced into the market [7j . We accompany our theoretical work with 
experimental evidence using a collinear setup for our type II periodically poled KTP 
(PPKTP) crystal. Our theory is well capable of reproducing the measurement results. 

The rest of the paper is organized as follows. In Sec. O we develop the theoretical 
model for our HSPS, and evaluate its second-order coherence properties as functions 
of source parameters in the ideal limit of infinitely high time resolutions. For this 
purpose, we use a heuristic continuous-mode analysis, whose validity is confirmed by an 
asymptotic discrete-mode analysis presented in the Appendix. Section [3] describes the 
experimental setup for our HSPS and the corresponding HBT interferometer, followed 
by our experimental results in Sec. [H There, we introduce our time-averaged coherence 
functions and their relation to the ideal figures. Section [5] concludes the paper. 



2. Theory of SPDC-based Single-Photon Sources 

The HSPS considered here consists of a parametric down-converter — driven by a cw 
pump at center frequency ujp producing cw signal (s) and idler (i) beams at center 
frequencies Ug and Ui = uUp — uUg^ respectively — followed by a single-photon detector 
on the idler beam. Here, we implicitly assume that signal and idler beams can be 
separated into two orthogonal spatial modes. In our experimental setup, this has 
been achieved by employing a type-II crystal, which creates signal and idler beams 
with orthogonal polarizations, along with a polarizing beam-splitter (PBS). Here, for 
simplicity, we suppress the spatial and polarization characteristics of signal and idler 
beams and represent them with scalar photon-units positive-frequency field operators 

m 

Ejit) = ^J du;Aj{u;)e-'-\ j = s,i, (1) 

where [Ej{t)^ Ej{u)] = 5{t — u) and Aj{uj) represents the corresponding output field 
operator in the frequency domain. It has been shown that, in the Heisenberg picture, 
the output field operators can be related to the vacuum-state field operators at the input 
to the crystal, A^-^{uj)^ j = s^i^ via the following Bogoliubov transformation j9j 

As{ujs +i^) = /i(^)irK + ^) + ^(^)^r^(^. - (2) 

A,{uj, -uj) = ^i{uJ)A'^{uJ^ -uj) + v{uj)A':\uJs + uj). (3) 
where — = 1. The joint state of signal and idler is a zero-mean Gaussian 

state whose only nonzero second-order moments are given by its temporal auto- and 
cross-correlation functions as follows P, 

{E]{t + T)E,{t)) = e-^-i?(r), i?(r) = J ^\u{u)\'e-, (4) 
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{E,{t + T)E,{t)) = (1 - 5,fc)e--''*--^-C7(T), C(r) = J '^v{uj)^i{uj)e-'-\ (5) 

where 5jk is the Kronecker delta function and j^k = s^i. 

In the low-gain regime, which is of interest to us, 1 and 

(^SPDc/5sPDc)sin^[^/(25sPDc)]/[^/(25sPDc)]^, where i?sPDC is the rate of photon 
generation for the signal/idler beam, and Sspdc is the bandwidth of the SPDC process 
pUj . In this regime, we have 

[ i?sPDc(l + rSsPDc) < r < 

R{t) = I i?sPDc(l - t5spdc) < t < ^ , (6) 

I elsewhere 

and 

I elsewhere ' 

where we assume that the difference in the speed of light for ordinary and extraordinary 
axes in the crystal has been compensated. In our experiment, B^pbc is on the order of 
THz, and i?sPDC is on the order of 1 MHz per milliwatt of pump power. 

In this paper, we calculate two coherence measures for our HSPS. The first 
figure quantifies the reliability of our heralding mechanism by looking at the temporal 
correlation between the signal and idler beams, and the second measure quantifies its 
capability to create one — and only one — photon per heralding event. In both cases, 
we first find the full temporal shapes of the above coherence functions in the ideal 
limit of infinitely high time resolutions and, then, later in Sec. [4] we will introduce our 
corresponding time-averaged parameters that we can measure in a typical experimental 
setup. Also, throughout the paper, we neglect the dark count effect and we assume that 
all employed photodetectors have unity quantum efficiencies. The latter assumption 
does not affect our measurement results because all correlation functions that we deal 
with in this paper have normalized forms. 

2.1. Signal-idler Temporal Correlation 

As a measure of temporal correlation between signal and idler, we obtain the degree of 
second-order coherence between the signal and the idler fields defined as follows, 

(2),, _ _ {Elit + T)E\{t)%{t)E,{t + r)) 



{Ei{t+T)E,{t^T)){E]mm 

|C(r)£_ (2). . ... 
= l + ^^=^.(r), (8) 

where, in the last step, we used the quantum form of the Gaussian moment-factoring 
theorem [9j by which we can reduce the fourth-order moment in the above equation to 
the sum of products of second-order moments, available from Eqs. (|4j) and ([5]), as follows 

P,,{t + r, t) = {El{t + T)Ej{t)E,{t)Es{t + r)) 
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= {Elit + r)Ej{t))mm{t + r)) 
+ {El{t + T)E,{t)){Ej{m{t + T)) 
+ {El{t + T)EAt + T)){Ej{m{t)) 

= R'{0) + \C{t)\'^P,,{t). (9) 

Here, Psi{t + r, t) is the coincidence rate for observing a signal photon at time t + t 
and an idler photon at time t, and from the above equation, it is only a function of r. 
In the low-gain regime, ^'^^^(O) ^ ^spdc/^spdc, which is inversely proportional to the 
probability of detecting a photon in a time interval of width At = 1/Bspbc- For our 
experimental setup, g^^\o) is on the order of 10^ and gfiij) has a narrow sub-picosecond 
width. These two properties witness ultrashort, highly correlated, twin wavepackets. 

2.2. Second-order Coherence Function for Heralded Signal 

The second coherence measure that we consider here is the degree of second-order 
coherence for the signal field, conditioned on observing an idler photocount at time 
t^, defined as follows 

(Bj(«l)B.(«l))p„{Bj((2)B.(«2))pm 

where (•)pm is the average over the post-measurement state assuming sufficiently high 
time resolution and unity quantum efficiency for the idler photodetector. 

To model the measurement on the idler field operator, we use a heuristic approach in 
which a photodetection event at time ti on the idler beam is modeled by the continuous- 
time measurement operator, [11], Ei{ti). In the Appendix, we employ a discrete-mode 
formalism for the same problem and show that in the asymptotic limit of infinitely 
high time resolution the results of the two methods converge. The post-measurement 
averaging, for any operator X, will then be given by 

(X)p„ = {^{U)X%{U))/{E\{U)%{U)). (11) 
The conditional coherence function in Eq. f fTOl ) can then be written as follows 

g^^\ti,t2\ti) = ^S\ti^h^ti)R{0) ^ ^^2^ 

Psi{ti^ ti)Psi{t2^ ti) 

where, using again the quantum version of the Gaussian moment-factoring theorem 
along with Eqs. (jlj) and ([5j), 

= i?(0) [R'iO) + \R{n,)\' + \Cin)\' + \C{r,)\'] 

+ 2^{C{n)C*{T2)R{Tu)} (13) 

is the multi-coincidence rate for finding signal photons at times ti and t2 and an idler 
photon at time t^. In the above equation, ri2 = ti — t2^ Ti = ti — ti^ and T2 = t2 — t^. 
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Finally, by plugging Eqs. f|T31) and into Eq. f|T21) . we find 

9fHtl,t2\ti) 



1 , 1 , |i^(ri2)lVi^^(0)-l 



^ 2^{(7(ri)(7-(r2)i?(ri2)} 

There are several interesting cases to be considered. First, let us look at the 
coherence function at the trigger time, i.e.. 

It is clear that if 5^5^^ (0) ^ 1 then gf\ti^ti\ti) as desired. In other words, the 
reliability of the heralding mechanism as well as the multiple-photon suppression are 
both guaranteed by the same condition i?^(0) <C |C(0)p. 

The second interesting case is when ti = ti but \t2\ = \t2 — ^ 2 At. In this case. 



+ 



23?{C(0)C*(r2)i?(-r2)} 



^ 1, (16) 

provided that g^^\o) ^ 1 and 9^^\t2) ~ 1- This implies that our HSPS has a coherence 
time on the order of At. 

Finally, let us consider the case when \ti = 3> 2At, i.e, when there is no 
correlation between the trigger time and the signal beam. In this case, 

gi'\t,MU) = 2 + ^^(1-^ 
9 si {^2) V ^^(^2), 

~ 2, (17) 

(2) 

provided that g^^ {T2) ^ 1, which prevails in the low-gain regime. This is in accord with 
the fact that in the SPDC process, in the lack of any triggering event, both signal and 
idler beams individually obey the thermal-state statistics, for which the second-order 
coherence function has a maximum value two 



3. Experimental Setup 

In this section, we describe the experimental setup used to demonstrate the theoretical 
findings from the previous section. Figure [1] presents the optical setup used for our 
HSPS along with the HBT interferometer used for the Qc measurement. A cw blue laser 
at center wavelength 405 nm pumps a type-II periodically-poled KTiOP04 (PPKTP) 
crystal. The crystal was from Raicol with a 10 |Lim period and its dimensions were 
1 X 2 X 10 mm^. The periodicity was chosen so that we would achieve creation of 
degenerate photon pairs at 810 nm slightly above the room temperature. The PPTKP 
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Figure 1. Experimental setup for our heralded single-photon source. A blue laser 
pumps a 1 X 2 X 10 mm^ PPKTP crystal to create signal and idler beams. The pump 
beam will be removed by using dichroic filters as well as interference filters. Signal and 
idler beams are split into different spatial modes by using a polarizing beam splitter 
(PBS). The idler beam is used as a trigger and the signal beam goes through a 50/50 
beam splitter for the ^^^^ measurement. 



crystal was inserted into a home-made oven made from copper and surrounded by PEEK 
plastic, which allowed us to reach the degeneracy at 39 °C, as shown in Fig. O with a 
stability of ±0.1 °C. Some optical elements were used to focus the laser, to reshape its 
spatial mode, and to collect the signal and idler beams. Dichroic filters were employed 
after the crystal to remove the pump beam. The crystal was cut for propagation along 
the X-axis in order to support mainly type-II SPDC for which the signal and idler photons 
have orthogonal polarizations. With the help of a PBS we could then deterministically 
split the two beams into two different spatial modes. A photodetection event on the 
idler beam heralds the presence of one or more photons on the signal beam, which goes 
through an HBT interferometer consisting of a 50/50 beam splitter, two interference 
filters, and two single-photon photodetectors. All photodetectors were single-photon 
counting modules from Perkin-Elmer with equal nominal quantum efficiencies of 0.4, 
dead-times of 45 ns, and time resolutions of 350 ps. Our interference filters had a 10 nm 
bandwidth centered at wavelength 810 nm. The measured photon count rate for our 
setup was up to around 850,000 counts/s in each channel at 50 mW pump power, with 
a signal-idler coincidence count rate amounting to at most about 10% of that value, 
which implies that the mode-matching was still far from perfect. Moreover, in order to 
avoid overloading the data acquisition hardware we occasionally chose to attenuate the 
down-converted beams with neutral density filters, which reduced the coincidence count 
rate even further without affecting the coherence properties of the source. 

The detection times for the signal and idler beams were recorded by a time-tagging 
card from Dotfast Consulting with a nominal temporal resolution of 156.25 ps. The 
time-tagging card streams the time tags to a computer by which we could calculate 
any single, double, or triple coincidence rates between the three channels (i, 5i, and 
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Figure 2. Tuning curve of the signal (Horizontal) and idler (Vertical) photon center 
wavelengths as a function of the crystal temperature. 
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Figure 3. Spectrum of the down-converted pair of photons at nominal degeneracy 
at 39 C with center-of-mass wavelengths of 809.82 nm (vertical) and 810.12 nm 
(horizontal). The different shapes and widths of the two spectra are to be expected 
from the material dispersion curve and the phase-matching relation, but may in part 
result from imperfect coupling to the multimode fibers used for spectroscopy and 
a nonuniform background during the measurement. For our purpose of building a 
heralded single-photon source, spectral indistinguishability is unimportant. 



52) in Fig. d] with a coincidence window that could be varied from 0.5 ns to 20 ns. The 
complete system of photodetectors, power supplies, time-stamping electronics and the 
USB interface fits in a 30 x 30 x 30 cm^ box. 

Historically, PPKTP crystals have mostly been used for type-I SPDC, i.e., identical 
polarization of the output photons, because the effective strength of the nonlinearity is 
lower for type-II than for type-I SPDC (for type-I, d^s = 10.7 pm/V and for type- 
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II, ds2 = 2.65 pm/V [12j). Nevertheless, as described previously, the signal and idler 
photons can be separated deterministically in the type-II case. For the second-harmonic 
generation (SHG), we found a conversion efficiency of about 0.03% at 80 mW pump 
power, which is close to the reported values in the literature. The main advantage of 
using PPKTP is that the poling enables collinear conversion via quasi-phase matching, 
which substantially improves the collection efficiency. 

Figure [3] shows a typical spectrum of the down-converted photons at 39° C. We 
used a 750 mm focal length spectrometer with 600 grooves/mm grating to obtain these 
spectra. We can clearly see the bimodal behavior of the twin photons as we select the 
spectrum in polarization. We can estimate a spectral bandwidth of AA = 5 and 7.5 nm 
corresponding to values previously reported in the literature [13j. 



4. Experimental results 

In this section, we report on our experimental results for the two coherence functions 
described in Sec. [2l In order to measure gf^ij)^ we approximate Psiir) in Eq. ^ by 
the rate of coincident events, Nsi{r)^ in which an idler photocount is observed at time 
t and a signal photocount is observed in the interval [t + r — Tcoin, t + r + Tcoin], where 
2T"coin is the width of our chosen coincidence window. Because of the photodetectors' 
time jitters, and neglecting dark counts throughout the paper, a photodetection event at 
time t only implies the existence of one or more photons in a neighborhood around time 
t. For simplicity, we assume that the detection time corresponding to a photon that hits 
the detector's surface at time t is uniformly distributed over the interval [t — Td^t + r^i], 
where Td is the time resolution of the photodetectors. We can then write the observed 
value for A^si(r) in terms of Psi(r) in the following way 

Nsiir) ~ / dr'Psiir'), 

where 



(18) 



Pdr) = j dUj dtsu{t,)u{ts - r)Ps^{ts - t,) (19) 

is the coincidence rate for detecting a signal photon(s) at time t+r and an idler photon(s) 
at time t, where u{t) = l/{2Td) if \t\ < r^i, and zero otherwise. 

Figure H] shows the experimental and the theoretical results for the time-averaged 
coherence function 

£\r)^Nsdr)/R\0) (20) 

for different values of pump power. Experimentally, i?^(0) was determined by the 
product of the signal and idler count rates. For the theoretical graphs, we used the 
low-gain correlation functions given by Eqs. ([6]) and ([71) with i?sPDC = 1.2 MHz per 
milliwatt of pump power and 5spdc = 3THz. It can be seen that 5^5^(0) drops as 
we increase the pump power, which is a direct result of multiple-photon contribution 
to the output. The peak value of ^'^^^(O) is also determined by the chosen coincidence 
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Figure 4. (Color online) Measurements (symbols) and theory predictions (lines) 
of the time-averaged coherence function gl^ (r) for the signal and idler photons at 
a chosen coincidence window of 0.78 ns. The low-gain regime theory curves are in 
striking agreement with the data using the following parameter values i^o/^pump ^ 
1.2-10^ pairs/ (s • mW), Td = 350 ps, and ^spdc = 3 THz. We measured the bandwidth 
by spectroscopy (see Fig. [3]). The pair production rate per pump power and the detector 
time resolution are approximated by subjective visual fitting to the above four data 
sets. 



window, here 0.78ns, because from Eqs. fl7l)-ffT9l). Nsi{r) has an almost fixed value for 
T G [— Tcoin + At + r^, Tcoin — At — T^] , inversely proportional to Tcoin- As we get farther 
from the center, the time-averaged coherence function drops to its minimum value one 
as expected. The theoretical graphs are in striking agreement with our experimental 
results, which clearly demonstrate the strong temporal correlation between signal and 
idler beams. 

To quantify multiple-photon suppression in our HSPS, we look at gi^\r) = 
gf\ti,ti + T\ti) = gi^\0,T\0). For an ideal HSPS, we expect that gi^\o) = 0. In 
our case, from Eq. ( 0^51 ). gi^\o) ^ 2 • 10~^ <C 1 at 15 MHz single-photon generation rate. 
However, again, we are only able to measure a time-averaged version of the coherence 
function by approximating Pgiir) with Nsi{T) as before and P^f (0,r,0) with A^i?(r), 
the count rate for a triple coincidence of an idler photodetection event at = 0, and 
two signal photodetection events at ti G [— Tcom, ^^coin] and t2 ^ [r — Tcom, ^ + ^^coin]- By 
accounting for the resolution of the three photodetectors involved in our measurement, 
we obtain 

NS\r) = - ^/ dh dt2PS\h,t2,0), (21) 



where 



(2Tcoin)^ 

pS\tut2,o) = I dti I dt,, I dt,Mti)u{ts, - h] 



X u{t,,-t2)PS\ts„t,„ti) (22) 



Characterizing heralded single-photon sources with imperfect measurement devices 11 




Figure 5. (Color online) The triple coincidence rate of having an idler photodetection 
event at time ti and signal photodetection events at times tgi and ts2 on, respectively, 
detectors si and s2 in Fig.[T]for a coincidence window of 2rcoin = 0.78 ns. Pump power 
for this measurement was 11.9 mW. 

is the multi-coincidence rate for detecting an idler photon(s) at time and two signal 
photons at times ti and t2- Figure [5] shows the experimental triple coincidence rate as 
a function of the two time differences. In this figure, the triple coincidence has been 
obtained by looking at the rate of an idler photodetection event at time ti and two 
signal photodetectoin events at times tgi and ts2 on, respectively, detectors Si and ^2 
in Fig. [H The peak at the center of the figure represents the contribution of multiple- 
photon pairs, and it is proportional to Pj^^^(0, 0, 0). The wall at tgi — = in Fig. [5] 
represents a coincidence event between the idler photon and one of the signal photons, 
detected by si, and it is proportional to ^^^^(0, r, 0), where r = ts2 — ti. Using Eq. f fT3l ). 
one can verify that the ratio between ^^^^(0, 0, 0) and Psf\o^ r, 0) is approximately given 
by (2i?2(0) + 4|C(0)n/(i?2(0) + |C(0)n ^ 4, where we assumed R\0) < \C{0)\^ and 
R{t) ^ C{r) ^ 0. The equivalent ratio obtained from Fig. [5] is about 2, which reflects 
the effect of time averaging in our analysis. 

Figure [6] shows our measurement results for the time- averaged conditional coherence 
function 

= NS\T)R{0)mmNdT)] (23) 

for different values of pump power, which result in different values for the observed 
central dip. Here, R{0) is obtained from the idler count rate in the experiment. The 
ringing structure in Fig. [6] is caused by double optical reflections [T4]. The graphs, 
nevertheless, exhibit the signature of a good SPS as the measured value of gc (0), at 
14 MHz single-photon generation rate, in Fig. [6l is below 0.03 for 2rcoin = 0.78 ns and 
Td = 0.35 ns. 
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Figure 6. (Color online) Measured (symbols) and calculated (lines) time- averaged 

(2) 

conditional coherence function gc \t). The theory lines were calculated using the 
same parameter values as in Fig. HI The purely statistical errors of our data are on 
the order of the symbol size in the figure and therefore not shown. As explained in a 
previous article [14J, photons that are reflected twice cause the apparent ringing. 
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Figure 7. (Color online) Measured (symbols) and calculated (lines) time- 
averaged conditional coherence function gi ' (0) versus pump power for three different 
coincidence windows. The theory lines are calculated using the same parameter values 
as in Fig. [H In the low-gain regime, there is a linear increase in g^c^ (0) versus pump 
power due to the multi-photon contribution to the down- converter output. 
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Figure 8. (Color online) Measured (symbols) and calculated (lines) time- averaged 
conditional coherence function gc (r) for three different coincidence windows. Here, 
the pump power was 11.9 mW corresponding to a pair production rate of approximately 
14 MHz. The theory lines are calculated using the same parameter values as in Fig. [H 
The theory curves only reproduce the data near the center of the dip and at very 
long delay times. This is to be expected because the finite time resolution and the 
shape of the spectrum were modeled with simplified rectangular and triangular shapes, 
respectively. 

By reducing the pump power we can reduce gi^"^ (0) almost arbitrarily at the expense 
of reducing the total count rate. This effect has been shown in Fig. [71 where we have 
plotted gP{0) versus the single-photon generation rate, i^sPDC, or equivalently, the 
pump power. In our experiment, each milliwatt of pump power corresponds to about 
1.2 million generated photon pairs per second. There is a linear growth in gi^\o) as a 
function of pump power, which exemplifies the contribution of multiple-photon states 
to the output in the low-gain regime. In this regime, from Eqs. f fMl ) and ([6])-([8]), 
gi^\0) ^ 2/^2^(0) ^ 2i?2(0)/|C(0)p = 2i?spDc/5sPDC, which is proportional to the 
pump power. The value of gi^\o) is also a function of coincidence window as shown 
next. 

In Fig. M one can see an example of how the conditional coherence function varies 
with the chosen coincidence window. Here, we have shown ^c^^(r) for three values of 
the coincidence window. It can be seen that the width of the central dip is almost given 
by 2rcoin- The measured value for gi^\o) goes down by choosing shorter coincidence 
windows. It does not, however, approach the actual value of gi^\o) so long as the 
detector time resolution Td ^ At. In order to make this point clearer, in Fig. M, 
we have plotted gi^\o) versus 2rcoin- It can be seen that, for Tcom <^ ^d^ 9c^\^) is 
determined by r^, whereas, for Tcoin ^ ^d) it is almost linearly increasing with Tcoin- Our 
theoretical treatment is again well capable of reproducing the measurement results. The 
graph shown in Fig. [9] exemplifies the fact that a single value for g^c\^) does not bear 
sufficient information to quantify the source performance. At a fixed rate, the interplay 
between the coincidence window and the time resolution of photodetectors must also 
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Figure 9. (Color online) Experimental (symbols) and theoretical (line) results for the 
minimum of the time-averaged conditional coherence function, ^c^^(O), as a function of 
the coincidence window 2rcoin using the same set of parameters as in Fig.|4]at a pump 
power of 11.9 mW. The dashed line is for ideal photodetectors {rd = 0). 

be accounted to give a proper figure of merit for an SPS. Eventually, the true value of 
gi^\o) can be obtained from Eq. f fTSl ) by estimating R{0) and C{0). This can be done 
by finding the parameters that can best reproduce all or a subset of graphs shown in 

Figs, m 

5. Conclusion 

In this paper, we theoretically and experimentally studied the coherence properties 
of heralded single-photon sources that use parametric down-conversion. We used the 
Gaussian characteristics of down-converted fields to analytically find the degree of 
second-order coherence between signal and idler fields as well as for the signal field, 
individually, when it is conditioned on the detection of an idler photon. Our theory 
is well capable of reproducing our experimental results, which demonstrated a high- 
quality source of sub-picosecond single photons. It also allowed us to study the impacts 
of the chosen coincidence window, the down- conversion parameters, and the resolution 
of photodetectors on the outcome. Such an analysis enables proper standardization of 
single-photon sources even with imperfect measurement devices. 
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Appendix, gi^'' calculation: discrete-mode formalism 



The analysis in this Appendix is based on looking at the system's behavior within a finite 
time interval or frequency band. In such cases, instead of working with the continuous- 
time field-operator formalism represented by Ej{t) and Aj{(jj)^ introduced in Sec. [21 we 
can deal with a discrete set of annihilation operators. Here, we first develop such a 
multi-mode but discrete representation for the field operators in time and frequency 
domains. We then use our new formalism to describe the system's initial state and the 
measurement on the idler beam, as well as to find the post-measurement state and the 
conditional coherence function. 

For a time interval of finite width T, such as [— T/2,T/2], the field operator in 
Eq. ([T]) can be written as [15j 

e'-^'m = J2hn'-^^^^^^^, te[-T/2,T/2]^ndj = s,i (24) 



where 



, . ,exp\2i7rnt/T] , , 

dte'-^'Ejit) , j = s,i. (25) 

-T/2 Vl 



Here, {exp[—27rint/T]/^/T}^ for integer n, forms an orthonormal set of basis functions 
that span all finite-energy functions over t G [—T/2, T/2]. The operator dj^n is the 
corresponding annihilation operator associated with the nth mode function, which 
represents a frequency band of eflFective width 27r5f = 27r/T around center frequency 
27m5f. Hence for T ^ At, the operators a^^^, satisfying [dj^n^al^] = Snm^jk^ for 
j^k G {s^i} and integers m and n, can describe the spectral behavior of the SPDC 
process with sufficient accuracy. 

Alternatively, one can span the spectral field operators in Eqs. ([2j) and ([3j) as follows 

+ uj,) = J2 5^. „!^^p2!^^1^, / ^ aj/27: G [-W, W] and j = s, t, (26) 

where 

? r .r-.^ . ,,exp[-2i7rn//(2iy)] 



j-w V2W 

Here, j = 5, i, is the annihilation operator associated with a time interval of effective 
width St = 1/{2W) centered at nSt. Again, if we choose W to be much larger than 
^SPDC, the operators bj^n can address the temporal behavior of the SPDC process with 
sufficient resolution. 

Assuming W > Bsfbc and TW ^ 1, the above dual pictures can be 
related to each other by plugging Eq. ([T]) into Eq. ( l25l ). and then, approximating 
J dfAj{Ljj + 27rf)exp{-2i7rft) by J^^ dfAj{u;j + 27rf) exp{-2i7rft). Then, with the 
help of Eq. ( l26l l and some algebra, one can obtain 



-WT ^ ^ 



Characterizing heralded single-photon sources with imperfect measurement devices 16 



- v-^ . exp\—2i7rnm/M] , , 

hn- 2^ Hm yfj , J = s,i (29) 

m=-WT 

where M = 2WT + 1, assumed to be integer, denotes the total number of modes 
considered for the description of the SPDC output. 

The above dual formalism enables us to analytically describe the initial state of 
the system, the measurement performed on the idler beam and the corresponding post- 
measurement state of the signal beam, as well as the coherence functions of our interest. 
The spectral representation given by {aj^n} allows us to describe \ipsi)^ the state of the 
system at the outcome of the parametric down-converter, explicitly in the following form 

\^s^) = 6^\i^s^,n). (30) 



where, from Eqs. ([2]) and ([3j), 

OO 

is the two-mode squeezed state associated with the joint state of harmonic oscillators 
represented by CLg^n and di-n- In the above equation, \k)aj^ is the fc-photon number 
state associated with a^^^, j = s^i^ and /x^ = ii{2Tm5f)^ Vn = iy{27m6f). 

Our calculations here mostly rely on an equivalent form of the above joint state, 
i.e., its Wigner characteristic function defined as follows 

= exp[-(|^„|2 - l/2)(|Cr + iC.r) + 25R{^„z.„CC}], (32) 

where D{dX) = ^^p[Cd^ — Co] is the displacement operator associated with the 
annihilation operator a. The main feature of the above characteristic function is its 
being Gaussian with respect to its complex arguments (g and Q. 

By using Eq. (I29l) . we can also find the Wigner characteristic function associated 
with the state of the temporal modes of the system. The one which is of interest to us 
for our future calculations is 

^J^fc A,^A,o j^^^^^^ 7s,^ 7i,o) = {D{hs^k, ls,k)D{hs^u^s,i)D{\Q, 7i,o)) 

exp[2i7r/cm/M] 



exp 



WT 

al 



m=—Wl ^ 



m=—WT ^ 
WT 

+ 7i,o Yl 4,J^-H.c. 

m=-WT 



= exp [-(i?o - l/2)(|7.,fcr + It.,/!' + iTi.oH 

+2^{t,nloCi}] , k I and \k-l\ < M, (33) 
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where H.c. denotes Hermitian conjugate, 

111 12 

Rn^ M ^ + R{nbt)im). n = -H^r...T^r,(34) 

WT 

Cr,= '^^^e-^™'^ ^C{n^^^ n = -WT---WT, (35) 

m=-WT 

and we have assumed that M is a sufficiently large prime number. The joint Wigner 
characteristic functions of any combination of signal and idler modes can similarly be 

calculated. In particular, xt^'"^''°(75,/c, 7i,o) = xt^'"^''''^''°(75,it, 0, 7^,0). 

The characteristic function in Eq. (I33l) has a Gaussian form and can tell us about 
the joint signal-idler state at different epochs of time. For instance, the joint state of bg^k 
and bi^o is entangled if and only if Ck ^ 0. That implies that, in the low-gain regime, 
modes represented by bg^k and bi^o are in separable states if and only if \k\ > At /6t. 
A click on the idler's photodetector at time ti = then has only correlation with 
photons appearing in the signal beam during [—At, At] interval. We clarify this issue 
by calculating the gi^\ti^t2\ti) below. 

Without loss of generality, we assume = 0, and, within our discrete-time 
formalism, we approximate gi^\ti^t2\0) by 

-7|1%¥t' (36) 

where the averaging is taken over the signal's post-measurement state, and k and / are 
integer numbers that satisfy ti G [{k — ^)St^ [k + \)5t) and ^2 ^ [{I — \)^t', + |)^^)- 

The post-measurement density operator, after a detection event on the idler mode 
represented by ^^^o, for our three modes of interest, represented by bg^k-, bs,h S^^o, for 
A; 7^ /, is given by [16j 

^det J 71 J 71 J 71 

MMKk. -ls^k)D{Ku -7s^i)D{ko. -7z,o)M;, (37) 

where f d'^a = d^a dQa, 
J J —00 J —00 ' 

Pdet = {Ai\M}Mi\ij,i) = 1 - l/i?o, (38) 

and the measurement operator Mj is defined as follows 

= 4,0- 10)6,06,0 (0|, (39) 

where lO)^^^ is the vacuum state and I^^^ is the identity operator associated with bi^o 
mode. The above measurement operator accounts for one or more idler photons in an 
interval of width 5t around ti = 0. 

Similarly, the post-measurement density operator for temporal modes represented 
by bs^k and 6^,0 is given by 

P^r^ = ^ [ ^xfc"'^"(7.,.,7.o)M,^(6.,,,-7^ (40) 

^det J 71 J 71 
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With the help of Eqs. (1371) and (HOl) and the Gaussian form of the characteristic 
functions, we can show that the coherence function in Eq. f l36l ) is given by: 



(41) 





\Rk 




\' + \Ck\' + \ 


\Ci\ 




{Q' + \ 


Ck\'){Q' + \Ci\^) 



^ 2i?o(i?o — ^)^{CkCiRk-i} 

2(i?o-i)|c,na|2 



(g2 + |Cfc|2)(g2 + |Q|2)' 

with Q"^ = Ro{Ro — 1)^ and k ^ I, and 



.S'(M10)^2-i|^. (42) 

One can verify that, in the hmit of Vl^ ^ oc, Eqs. ( HB and (Il2l) converge to Eq. (fT2l) . 
This is because, in this regime, 5t ^ 0, hence the discrete-time annihilation operators 
approach to the continuous-time field operators. This proves that the heuristic approach 
that we employed in the previous section is indeed valid and gives us the correct result 
if the idler's photodetector has zero time jitter. One can also verify that the above 
equations reproduce all special cases we considered previously. 
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